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A measurement of the ZZ production cross section in proton-proton collisions at -yi = 7 TeV 
using data recorded by the ATLAS experiment at the LHC is presented. In a data sample cor- 
responding to an integrated luminosity of 1.02 fb~^ collected in 2011, 12 events containing two 
Z boson candidates decaying to electrons and/or muons are observed. The expected background 
contribution is 0.3 ± 0.3(stat.)^g'3(syst.) events. The cross section measured in a phase-space re- 
gion with good detector acceptance and for dilepton masses within the range 66 GeV to 116 GeV is 
'^^'z^i+e- e+e- ~ 19.4j|g'2(stat.)]!|g'7(syst.) ± 0.7(lumi.) fb. This result is then used to derive the to- 
tal cross section for on-shell ZZ production, a^zz = 8.5l2'3(stat.)^g'3(syst.)±0.3 (lumi.) pb, which 
is consistent with the Standard Model expectation of 6.51q'2 pb calculated at the next-to-leading 
order in QCD. Limits on anomalous neutral triple gauge boson couplings are derived. 

PACS numbers: 14.70.Hp, 12.15.Ji, 12.60. Cn, 13.85. Qk 



The production of pairs of Z bosons at the LHC is of 
great interest since it provides an excellent opportunity 
to test the predictions of the electroweak sector of the 
Standard Model at the TeV energy scale; moreover it is 
the irreducible background to the search for the Higgs 
boson in the H — ZZ decay channel. In the Standard 
Model, ZZ production proceeds at leading order (LO) 
via t-channel quark-antiquark interactions; the ZZZ and 
ZZj neutral triple gauge boson couplings (nTGCs) are 
absent, hence there is no contribution from s-channel qq 
annihilation at tree level. At the one-loop level, fermion 
triangles generate nTGCs of 0(10"'*) [T]. Many mod- 
els of physics beyond the Standard Model predict values 
of nTGCs at the level of lO--* to lO'^ gj. The signa- 
ture of non-zero nTGCs is an increase of the ZZ cross 
section at high ZZ invariant mass and high transverse 
momentum of the Z bosons [3] . ZZ production has been 
studied in e+e" colhsions at LEP [U [5] and in pp col- 
lisions at the Tevatron [BJ [7] . No deviation of the mea- 
sured cross section from the Standard Model expectation 
has been observed, and limits on anomalous nTGCs have 
been set 5, 6J. 

This letter presents the first measurement of ZZ [8 
production in proton-proton collisions at a centre-of-mass 
energy ^/s of 7 TeV, and limits on the anomalous nTGCs. 
The cross section for on-shell ZZ production (i.e. in 
the zero-width approximation) is predicted at next-to- 
leading order (NLO) in QCD to be 6.5toi pb 0; this 
includes a ~6% contribution from gluon fusion. Candi- 
date ZZ events are reconstructed in the ZZ — >■ £~^£~£~^£^ 
decay channel, where £ can be an electron or muon. Al- 
though this channel constitutes only ~0.5% of the total 
ZZ cross section, its final state with four high transverse- 
momentum, isolated leptons has a very high expected 
signal to background ratio of ^30. 

To reduce systematic uncertainties, the cross section is 
measured within a phase-space that corresponds closely 
to the experimental acceptance; this is termed the 'fidu- 
cial' cross section. The fiducial phase-space definition 



requires the invariant mass of both lepton pairs to be be- 
tween 66 GeV and 116 GeV and all four leptons to be 
within the pseudorapidity |10| range \ri\ < 2.5 and have 
transverse momentum px > 15 GeV. The four-momenta 
of all photons present after the simulation of the parton 
shower which are within AR = \/ Acj)^ + Arf < 0.1 of a 
lepton are summed into the four momentum of that lep- 
ton. The total ZZ cross section in the on-shell approx- 
imation is obtained from the fiducial cross section using 
the known Z — ^ £^ £~ branching ratio and a correction 
factor for the kinematic and geometrical acceptance. 

Anomalous nTGCs for on-shell ZZ production can be 
parameterized by two CP-violating {fX) and two CP- 
conserving {fY) complex parameters {V = Z,j) which 
are zero in the Standard Model [3|. To ensure partial- 
wave unitarity, a form-factor parameterization is intro- 
duced to cause the couplings to vanish at high parton 
centre-of-mass energy fi" = /j^/(l-fs/A^)". Here, A 
is the energy scale at which physics beyond the Standard 
Model will be directly observable, are the low-energy 
approximations of the couplings, and n is the form-factor 
power. Following Ref. [3], n = 3 and A = 2 TeV are 
chosen, so that expected limits are within the values pro- 
vided by unitarity at LHC energies. The results with 
energy cutoff A = oo are also presented as a comparison 
in the unitarity violation scheme. 

The ATLAS detector [TT] consists of inner tracking de- 
vices surrounded by a superconducting solenoid, electro- 
magnetic and hadronic calorimeters and a muon spec- 
trometer with a toroidal magnetic field. The inner detec- 
tor, in combination with the 2 T field from the solenoid, 
provides precision tracking of charged particles for |7/| < 
2.5. It consists of a silicon pixel detector, a silicon strip 
detector and a straw tube tracker that also provides tran- 
sition radiation measurements for electron identification. 
The calorimeter system covers the pseudorapidity range 
I77I < 4.9. It is composed of sampling calorimeters with 
either liquid argon (LAr) or scintillating tiles as the active 
media. In the region |?7| < 2.5 the electromagnetic LAr 



2 



calorimeter is finely segmented and plays an important 
role in electron identification. The muon spectrometer 
has separate trigger and high-precision tracking cham- 
bers which provide muon identification and measurement 
in |?7| < 2.7. 

A three-level trigger system selects events to be 
recorded for offiine analysis. The events used in this 
analysis were selected with single-lepton triggers with 
nominal transverse momentum thresholds of 20 GeV for 
electrons and 18 GeV for muons. The efficiencies of the 
single-lepton triggers have been determined as a function 
of lepton pt using large samples oi Z ^ £^£^ events. 
The trigger efficiency for events passing the offiine se- 
lection described below is 99.9% with an uncertainty of 
0.1%. 

This measurement uses a data sample of proton-proton 
collisions at a/s = 7 TeV recorded between February and 
June 2011. After data quality requirements, the total 
integrated luminosity used in the analysis is 1.02 fb~^. 
The integrated luminosity uncertainty is 3.7% |12| . 

Events are required to contain a primary vertex formed 
from at least three associated tracks. The vertex with 
the largest sum of the computed from the associated 
tracks is selected as the primary vertex. 

Signal events are characterized by four high-pT, iso- 
lated electrons or muons, in three channels: e~^e~e^e~, 



and e+e fj,~^^ 



Lepton candidates are re- 
quired to be consistent with originating from the primary 
vertex. Muons are identified by matching tracks (or track 
segments) reconstructed in the muon spectrometer to 
tracks reconstructed in the inner detector [TH] . Their mo- 
mentum is calculated by combining the information from 
the two systems and correcting for the energy deposited 
in the calorimeters. Only muons with > 15 GeV and 
1 77 1 < 2.5 are considered. In order to reject muons from 
the decay of heavy quarks, isolated muons are selected 
by requiring the scalar sum of the transverse momenta 
(Spt) of other tracks with px > 1 GeV inside a cone of 
size AR = 0.2 around the muon to be no more than 15% 
of the muon px- The overall reconstruction, identifica- 
tion and isolation efficiency, measured in data using a 
large sample of Z A^^M^ events, varies as a function 
of PT from 92% at 15 GeV to 95% at 45 GeV. 

Electrons are reconstructed from a cluster in the elec- 
tromagnetic calorimeter matched to a track in the in- 
ner detector [13]. Electron candidates are required to 
pass the 'medium' identification criteria described in 
Ref. [13 , to have a transverse momentum (measured in 
the calorimeter) of at least 15 GeV and a pseudorapidity 
of |?7| < 2.47. They must be isolated, using the same 
criterion as for muons, calculating the Epx around the 
electron track. Electron candidates within AR = 0.1 of 
any selected muon are rejected, and if two electron candi- 
dates are within AR = 0.1 of each other the one with the 
lower px is rejected. The overall reconstruction, identifi- 
cation and isolation efficiency varies as a function of px 
from 63% at 15 GeV to 81% at 45 GeV. 

Selected events are required to have exactly four lep- 



tons, and to have passed a single- muon or single-electron 
trigger. To ensure high trigger efficiency, at least one 
of these leptons must have px > 20 GeV (25 GeV) for a 
muon (electron) and match to a muon (electron) recon- 
structed online by the trigger system within AR < 0.1 
(0.15). 

Same-fiavour, oppositely-charged lepton pairs are com- 
bined to form Z candidates. An event must contain 
two such pairs. In the e^e^e+e~ and fi'^ fx^ fi^ fi^ chan- 
nels, ambiguities are resolved by choosing the pairing 
which results in the smaller value of the sum of the two 
\m£+i~ — mz\ values. Figure [l] shows the correlation be- 
tween the invariant mass of the leading (higher px) and 
the subleading (lower px) lepton pair. The events cluster 
in the region where both masses are around mz- Events 
are required to contain two Z candidates with invariant 
masses satisfying 66 GeV < mi+i- < 116 GeV. 
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FIG. 1. The mass of the leading lepton pair versus the mass 
of the subleading lepton pair. The events observed in the 
data are shown as solid circles and the ZZ signal prediction 
from simulation as boxes. The large dashed box indicates the 
signal region defined by the requirements on the lepton-pair 
masses. 



The reconstruction efficiency for ZZ events is deter- 
mined from a detailed Monte Carlo simulation. The LO 
generator Pythia [14^ with the MRST modified LO par- 
ton density function (PDF) set [TS] is used to model 
pp — > ZZ — l'^t~t^t~ events, where £ includes elec- 
trons, muons and r leptons. The Pythia simulation 
includes the interference terms between the Z and 7* 
diagrams; the mass threshold for the Z/7 boson is set to 
12 GeV. The detector response is simulated [16] with a 
program based on GEANT4 [TT . Additional inelastic pp 
events are included in the simulation, distributed so as 
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to reproduce the number of collisions per bunch-crossing 
in the data. The simulation is also corrected with scale 
factors, and the lepton momentum resolution adjusted, 
to reproduce the lepton reconstruction and identification 
efficiencies measured in data. 

The overall efficiencies of the reconstruction and se- 
lection criteria for events generated within the fidu- 
cial phase-space are (40±3)%, (79±2)% and (57±2)% 
for e+e^e+e^, and e'^ e~ fx'^ fi~ respectively. 

The dominant systematic uncertainties arise from elec- 
tron identification (6.6% in the e+e^e+e^ final state, 
3.1% in the e^e~ pL^^~ final state) and from the muon 
reconstruction efficiency (2.0% in ^'^ ii~ jX^ ^~ and 1.0% 
in ji^ 11^). 

Background to the ZZ signal originates from events 
with a Z (or W^) boson decaying to leptons plus addi- 
tional jets or photons {W/Z -f X), from top-quark pro- 
duction and from other diboson final states. Such events 
may contain electrons or muons from the decay of heavy- 
flavoured hadrons, or muons from in-flight decay of pi- 
ous and kaons; jets or photons may be misidentified as 
electrons. The majority of these background leptons are 
rejected by the isolation requirement. 

To estimate the background contribution from four- 
lepton events in which one lepton originates from a jet, a 
sample of events containing three leptons passing all se- 
lection criteria plus one 'lepton-like jet' is identified; such 
events are denoted Mj. For muons, the lepton-like jets 
are muon candidates that fail the isolation requirement. 
For electrons, the lepton-like jets are clusters in the elec- 
tromagnetic calorimeter matched to inner detector tracks 
that fail either or both of the full electron selection and 
the isolation requirement. The events are otherwise re- 
quired to pass the full event selection, treating the lepton- 
like jet as if it were a fully identified lepton. This event 
sample is dominated hy Z + X events. The background is 
then estimated by scaling this control sample by a mea- 
sured factor / which is the ratio of the probability for 
a jet to satisfy the full lepton criteria to the probability 
to satisfy the lepton-like jet criteria. The background in 
which two selected leptons originate from jets is treated 
similarly, by identifying a data sample with two leptons 
and two lepton-like jets; such events are denoted Ujj. To 
avoid double counting in the background estimate, and 
to take into account the expected ZZ contribution in the 
control region, N{ZZ), the total number of background 
events iV(BG) is calculated as: 

iV(BG) = N{mj) X / - N{e£jj) xf- N{ZZ). (1) 

The factor / is measured in a sample of data selected 
with single-lepton triggers with criteria applied to sup- 
press isolated leptons from and Z bosons, and cor- 
rected for the remaining small contribution of true lep- 
tons using simulation. It is measured independently in 
T] and pt and the values combined assuming they are 
uncorrelated. A similar analysis is performed on Monte 
Carlo simulation of background processes; the larger of 
the statistical uncertainty on / determined from the data 



TABLE I. Summary of observed events in the data, total 
background contributions and expected signal in the individ- 
ual four-lepton and combined channels. The quoted uncer- 
tainties represent 68.3% confidence intervals; the first is sta- 
tistical while the second is systematic. The uncertainties on 
the integrated luminosity (3.7%) and the theoretical ZZ cross 
section are not included. 
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and the difference between data and simulation is taken 
as the systematic uncertainty in each px (or 77) bin. This 
results in a systematic uncertainty which varies as a func- 
tion of pt from 57% (85%) at 15 GeV to 55% (77%) at 
45 GeV for electrons (muons). 

The numbers of expected and observed events after 
applying all selection criteria are shown in Table |lj The 
expected number of signal events is determined from the 
Pythia simulation normalized to the NLO calculation 
using MCFM m with the MSTW2008 [18] NLO PDF set. 
The normalization factor, calculated within the phase- 
space of the fiducial cross section measurement, is 1.41. 
The expected numbers of signal events include contribu- 
tions of 1.6% from ZZ — > £^£^£^£^ events generated 
outside the fiducial phase-space and 0.3% from events 
where one of the Z bosons decays to r leptons. Twelve 
ZZ candidates are observed in data, with a background 
expectation of 0.3 ± 0.3(stat.)^Q'3(syst.), corresponding 
to a p-value of 10^^ equivalent to a one-sided Gaussian 
significance of 5a. In the four- muon channel 8 events are 
observed where 3.3^0 3 signal plus background events are 
expected. The probability of the expected number fluc- 
tuating up to 8 or more is 3.2%. 

The transverse momentum distribution and the invari- 
ant mass distribution of the combined four-lepton system 
for the selected candidates are shown in Fig. [2] 

The ZZ fiducial cross section is determined using a 
maximum likelihood fitting method to combine the three 
four-lepton channels. The systematic uncertainties are 
included in the fitting procedure as nuisance parameters. 
The measured fiducial cross section is: 

^I'^z^i+e-i+i- = 19.4i|^(stat.)+°;^(syst.)±0.7(lumi.)fb, 

where £^£^£^£^ refers to the sum of the e"'"e~e+e~, 
e+e^/i+/i^ and 11^ fi^ fi'^ i-i^ final states. The total cross 
section is determined similarly, correcting for the known 
Z — >■ £^£~ branching ratios and the acceptance of the 
fiducial phase-space. This acceptance, calculated at NLO 
using MCFM version 6.0 with the MSTW2008 PDF set, 
is 0.507 ± 0.009, where the error arises primarily from 
PDF uncertainties with a 1% contribution from QED ra- 
diative corrections and off-shell Z/j* effects evaluated 
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FIG. 2. (a) Transverse momentum ^ and (b) invariant mass Mzz of the four-lepton system for the selected events. The 
points represent the observed data and the histograms show the signal prediction from simulation. The shaded band on 
each histogram shows the combined statistical and systematic uncertainty on the signal prediction. The predicted number of 
background events from the data-driven background estimate is indicated on the plot. 



TABLE II. One dimensional 95% confidence intervals for 
anomalous neutral gauge boson couplings, where the limit 
for each coupling assumes the other couplings fixed at their 
Standard Model value. Limits are presented for form factor 
scales of A = 2 TeV and A = cxd and include both statistical 
and systematic uncertainties; the statistical uncertainties are 
dominant. 
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oo 


[-0.08,0.08] 


[-0.07,0.07] 


[-0.08,0.08] 


[-0.07,0.07] 



from POWHEG BOX ^9j. The measured value of the 
total on-shell ZZ cross section is: 

af^z = 8.5^2^ (stat.) (syst.) ±0.3 (lumi.) pb. 

The result is consistent within errors with the NLO 
Standard Model total cross section for this process of 
6.5t°i pb M- 

Limits on anomalous nTGCs are determined using the 
total number of observed events only. The ZZ produc- 
tion yield dependency on couplings is parameterized us- 
ing fully simulated events generated with SHERPA [20] 
subsequently reweighted using the leading-order matrix 
element [3] within the framework of Ref. [5T]. The 
reweighting procedure uses simulated samples with Stan- 
dard Model as well as non-Standard Model coupling val- 
ues to ensure adequate coverage of all kinematic regions. 
One dimensional 95% confidence intervals for the anoma- 
lous nTGCs are determined using a maximum profile 



likelihood fit to the observed number of events. The sys- 
tematic errors are included as nuisance parameters. The 
resulting limits for each coupling, determined assuming 
real couplings and with the other couplings fixed at their 
Standard Model value, are listed in Table [HI The present 
results are dominated by statistical uncertainties: limits 
derived using statistical uncertainties alone differ from 
those in Table [TTl by less than 0.01. These limits are com- 
parable with, or are more stringent than, those derived 
from measurements at LEP |5j and the Tevatron [6]; it 
should be noted that limits from LEP do not use a form 
factor, and those from the Tevatron use A = 1.2 TeV. 

In summary, the ZZ production cross section has been 
measured in proton-proton collisions at y/s = 7 TeV us- 
ing the ATLAS detector. Both the fiducial cross section 
within the detector acceptance and the total cross section 
have been determined. The latter is in agreement with 
the Standard Model expectation. Limits on anomalous 
nTGCs have been derived. 
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